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NOMENCLATURE 



English Symbols 



A 


Area (ft^) 


A dr 


2 

Drum Downcomer Total Cross Sectional Area (ft ) 


A ds 


2 

Header Downcomer Total Cross Sectional Area (ft ) 


AMB 


2 

Main Bank Heat Transfer Surface Area (ft ) 


A 

r 


2 

- Main Bank Riser Total Cross Sectional Area (ft ) 


A 

s 


2 

- Screen Riser Total Cross Sectional Area (ft ) 


Be 


Becker's Factor 


C 


A General Constant 


C £ 


Specific Heat of Liquid (Btu/lbm/°F) 


C m 


Specific Heat of Tube Metal (Btu/lbm/°F) 


C P 


Specific Heat (Btu/lbm/°F) 


C R1 


Specific Heat of Flue Gas in the Furnace 
(Btu/lbm/°F) 


C Rla 


Specific Heat of Flue Gas at the Furnace 
Exit (Btu/lbm/°F) 


C R2 


Specific Heat of Flue Gas at Superheater Exit 
(Btu/lbm/°F) 


C R3 


Specific Heat of Flue Gas at Main Bank Exit 
(Btu/lbm./°F) 


C st 


Constant used in Rohsenow Correlation 


C 3- C ll 


Constants Used in Various Equations 


d 


Diameter (ft) 


D, 

dr 


- Sum of the Diameters of the Drum Downcomers 
(ft) 


D, 

ds 


- Sum of the Diameters of the Screen Downcomers 
(ft) 
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D - Sum of the Diameters of the Main Bank Tubes 

r (ft) 

D - Sum of the Diameters of the Screen Tubes (ft) 

s 

EB - Total Heat Input to the Furnace (Btu/sec) 

EDI - Heat Transferred from the Desuperheated 

Steam to the Tubes (Btu/sec) 

ED2 - Heat Transferred from the Superheater Tubes 

to the Steam (Btu/sec) 

ER1 - Heat Transfer Rate from the Flue Gas to the 

Screen Tubes (Btu/sec) 

ER2 - Heat Transfer Rate from the Flue Gas to the 

Superheater (Btu/sec) 

ER3 - Heat Transfer Rate from the Flue Gas to the 

Main Bank (Btu/sec) 

ER4 - Heat Transfer Rate from the Flue Gas to the 

Economizer (Btu/sec) 

E - Heat Transfer Rate to the Riser Liquid 

wr (Btu/sec) 

EWl - Heat Transfer Rate from the Economizer Tubes 

to the Feed Water (Btu/sec) 

EW2 - Heat Transfer Rate from the Desuperheater 

Tubes to the Drum (Btu/sec) 

EW3 - Heat Transfer Rate to the Screen Riser Liquid 

v (Btu/sec) 

EW4 - Heat Transfer Rate to the Main Bank Riser 

Liquid (Btu/sec) 

f - Friction Factor 

F, - Drum Downcomer Friction Factor 

dr 

F , - Heater Downcomer Friction Factor 

ds 

F - Main Bank Riser Friction Factor 

r 

F - Screen Riser Friction Factor 

s 

2 

g - Acceleration of Gravity - 32.2 ft/sec 

G - Mass Flow Rate (lbm/sec) 
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GB 


- 


Fuel Flow Rate (lbm/sec) 


GDI 


- 


Total Steam Flow Rate from Boiler (lbm/sec) 


GD2 


- 


Desuperheated Steam Flow Rate (lbm/sec) 


GD3 


- 


Steam Condensing from Vapor to Liquid in the 
Steam Drum (lbm/sec) 


GD4 


- 


Superheated Steam Flow (lbm/sec) 


GL 


- 


Air Mass Flow Rate (lbm/sec) 


GR 


- 


Total Flue Gas Flow Rate (lbm/sec) 


GS 


- 


Screen Riser Circulation Flow Rate (lbm/sec) 


Gr 


- 


Riser Circulation Flow Rate (lbm/sec) 


Gw 


- 


Downcomer Circulation Flow Rate (lbm/sec) 


G0 




Main Bank Riser Mass Flow Rate (lbm/sec) 


Gwr 


- 


Drum Downcomer Circulation Flow Rate 
( lbm/sec) 


Gws 


- 


Header Downcomer Circulation Flow Rate 
( lbm/sec) 


HDl 


- 


Enthalpy of Saturated Steam (Btu/lbm) 


Hfg 


- 


Heat of Vaporization (Btu/lbm) 


HO 


- 


Heating Value of Fuel Oil (Btu/lbm) 


HNB 


- 


Non-boiling Height (ft) 


HNBR 


- 


Main Bank Riser Non-boiling Height (ft) 


HNBS 


- 


Screen' Riser Non-boiling height (ft) 


HW1 


- 


Enthalpy of Saturated Water (Btu/lbm) 


HW2 


- 


Enthalpy of Downcomer Liquid (Btu/lbm) 


HW3 


- 


Enthalpy of Feed at Economizer Inlet 
(Btu/lbm) 


HW4 


- 


Enthalpy of Feed at Economizer Outlet 
(Btu/lbm) 


HW6 


- 


Enthalpy of Water in Lower Drum (Btu/lbm) 
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Hxr 


Effective Height of Main Bank/Drum 
Downcomer Circuit (ft) 


Hxs 


Effective Height of Screen/Header Downcomer 
Circuit (ft) 


k 


Thermal Conductivity (Btu/hr-ft-°F) 


K 


a general constant 


Ks 


Equivalent Sand Roughness (in) 


KR1-KR4 - 


Constants Used in the Gas Heat Transfer Equations 


l 


Length (ft) 


L 


Tube Length (ft) 


LB 


Boiling Length (ft) 


Ldr 


Total Length of Drum Downcomer Tubes (ft) 


Lds 


Total Length of Header Downcomer Tubes (ft) 


Ldlr 


Avg. Length of Drum Downcomer Tubes (ft) 


Ldls 


Avg. Length of Header Downcomer Tubes (ft) 


Lr 


Total Length of Main Bank Tubes (ft) 


Lrl 


Avg. Length of Main Bank Tubes (ft) 


Ls 


Total Length of Screen Tubes (ft) 


Lsl 


Avg. Length of Screen Tubes (ft) 


LMTD 


Log Mean Temperaturfe Difference (°F) 


M 


Mass (lbm) 


MM1 


Mass of Screen Tube Metal (lbm) 


MM2 


Mass of Superheater Tube Metal (lbm) 


MM3 


Mass of Riser Tube Metal (lbm) 


MM4 


Mass of Economizer Tube Metal (lbm) 


MM5 


Mass of Desuperheater Tube Metal (lbm) 


MRl 


Mass of Flue Gas in Furnace (lbm) 


MWl 


Mass of Water in Steam Drum (lbm) 
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p 

PD 

PDl 

PD2 

PD2 

Pr 

Psat 

PW1 

PW2 

q 

R 

Ral 

Ra2 

RD1 

Re 

RW 

RS 

R0 

t 

T 

Tsat 

TD 

TD1 

TD2 



Pressure (psia) 

Pressure of Steam Volume in Steam Drum (psi) 

Saturation Pressure Corresponding to Enthalpy 
of Water in the Steam Drum (psi) 

Pressure at Superheater Outlet (psia) 

Pressure at Desuperheater Outlet (psia) 

Prandtl Number 

- Saturation Pressure (psia) 

Pressure in Water Drum (psi) 

- Pressure in Screen Header (psi) 

Heat Transfer Rate (Btu/sec) 

Tube Radius (in) 

Average Density in the Main Bank Risers 
(lbm/ft 3 ) 

Average Density in the Screen Risers 
(lbm/ft 3 ) 

Density of Saturated Steam (lbm/ft 3 ) 

Reynolds Number 

3 

Density of Saturated Water (lbm/ft ) 

- Density of Two-Phase Mixture at Screen 
Outlet (lbm/ft 3 ) 

Density of Two^Phase Mixture at Main Bank 
Outlet (lbm/ft ) 



Time (sec) 

- Temperature (°F) 

- Saturation Temperature (°F) 

Temperature of Steam Drum Vapor (°F) 

Temperature of Saturated Steam Corresponding 
to Drum Enthalpy (°F) 

Superheater Outlet Temperature (°F) 
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TD3 

TMl 

TM2 

TM3 

TM4 

TM5 

TR1 

TR2 

TR3 

TR4 

TR5 

TRa2 

TRa3 

TRa4 

TW3 

TW4 

To 

U b 

U m 

V 

VB 

VDl 

VENT 

VNDR 

VR 

VS 



Desuperheater Outlet Temperature (°F) 

Screen Metal Temperature (°F) 

Superheater Metal Temperature (°F) 

Main Bank Metal Temperature (°F) 

Economizer Metal Temperature (°F) 

Desuperheater Metal Temperature (°F) 

Temperature of Flue Gas in Furnace (°F) 

Temperature of Flue Gas Leaving Furnace (°F) 

Temperature of Flue Gas Leaving Superheater (°F) 

Temperature of Flue Gas Leaving Main Bank (°F) 

Temperature of Flue Gas Leaving Economizer (°F) 

Average Flue Gas Temperature Across the 
Superheater (°F) 

Average Flue Gas Temperature Across the 
Main Bank (°F) 

Average Flue Gas Temperature Across the 
Economizer (°F) 

Feed Inlet Temperature (°F) 

Feed Outlet Temperature (°F) 

Ambient Temperature 

Bubble Rise Velocity (ft/sec) 

Mean Velocity (ft/sec) 

Volume (ft^) 

3 

Boiling Volume (ft ) 

3 

Volume of Steam in Steam Drum (ft ) 

Percent Throttle Valve Opening 

3 

Water Drum Volume (ft ) 

3 

Main Riser Bank Fluid Volume (ft ) 

3 

Screen Riser Bank Fluid Volume (ft ) 
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V 



V 


Specific Volume (ft'Vlbm) 


X 


Quality (percent) 


X0 


Quality at Main Bank Riser Outlet (percent) 


XS 


Quality at Screen Riser Outlet (percent) 


Xe 


Quality at Riser Exit (percent) 


Z1 


Drum Downcomer Bend Loss Factor 


Z2 


Drum Downcomer Exit Loss Factor 


Z3 


Main Bank Riser Bend Loss Factor 


Z4 


Main Bank Riser Exit Loss Factor 


Z5 


Screen Riser Bend Loss Factor 


Z6 


Screen Riser Exit Loss Factor 


Z7 


Header Downcomer Bend Loss Factor 


Z8 


Header Downcomer Exit Loss Factor 



Greek Letters 



V 


Viscosity (lbm/ft-sec) 


p 


Density (lbm/ft^) 


a - 


Surface Tension (lbf/ft) 


p Dal 


Average Density in Superheater (lbm/ft^) 


p Da2 


Average Density in Desuperheater (lbrn/ft"^) 



Subscripts 



a 


Average 


iff 


Liquid 


g 


Vapor 


M,m - 


Metal 
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r 

s 

dr 

ds 

D 

R 

W 



Far Field Average 
Main Bank 
Screen 

Drum Downcomer 
Header Downcomer 
Steam 

- Flue Gas 
Water 
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I. 



INTRODUCTION 



The increasing attractiveness and efficiency of computer- 
aided design techniques have established a need for 
reasonably accurate computer models of all types of 
engineering systems. The design of multivariable control- 
lers, in particular, requires the availability of state 
information about a system which is not available from 
conventional transfer function simulations. A number of 
computer models of marine boilers have been proposed 
[1,2, 3, 4]. All of these efforts share one or more of the 
following features: 

a. Lumped parameters are used. 

b. Bubble formation effects with respect to shrink and 
swell are not modeled. 

c. The furnace screen and main riser bank are lumped 
together as are their respective downcomers [1,3,4]. 

d. The model is linearized using a Taylor expansion 
or perturbation technique [3] . 

Use of lumped parameters is dictated by the need to maintain 
computational efficiency, particularly when the model is 
used in control system design and optimization schemes. 

In the past, effects due to vapor formation in the drum 
have not been included due to the lack of a suitable model. 

A simple model based on the Rohsenow correlation [5] was 
developed for use in the current simulation. However, its 
contribution to shrink and swell effects was negligible and 
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it was eliminated from the equation set. Regarding the 
tube banks, Haraguchi and Akasara [2] split out the various 
riser and downcomer loops, including intermediate headers. 
Considering the general level of uncertainty inherent in 
the correlations used in this model, this appears to be an 
excessive degree of complication. A compromise position 
was taken in the present effort in that the screen and main 
riser bank and their associated downcomers were modeled 
separately, but intermediate headers within the screen 
circuit were not addressed individually since their configura- 
tion varies widely from one boiler to another. 

The objective of this thesis was to develop a general, 
non-linear model of a Navy boiler for use in simulations 
covering the full operating range of the propulsion plant. 
However, since a linearized version is normally needed for 
multivariable controller design, CENTRL, a computer code 
developed locally for use with CSMP models, could be used 
to generate linearized state space representations described 
by the matrix equations: 



This approach would produce linearized representations 
covering the entire boiler operating range in a piecewise 
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fashion, rather than attempting to extend the range of 
applicability of a single state space or linear model. 

Since this author found the notation used in Reference 
[1] particularly readable, it was adopted for use in this 
document. 
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II. DESIGN CONSIDERATIONS 



A. D-TYPE BOILER 

The single furnace (D-Type) , uncontrolled-superheat 
boiler was chosen as the basis for the model since it is 
the predominant design in the current U.S. Navy inventory. 
The nominal operating pressures range from 600 to 1200 psi, 
and all designs are conceptually similar. A basic schematic 
is shown in Figure 1. These boilers are typified by 
relatively fast response times for use in high speed 
maneuvering vessels. This puts heavy emphasis on properly 
designed control systems. Manual control of these boilers 
is essentially infeasible except under the most favorable 
circumstances, e.g. steady steaming at low speeds. Current 
conventional control systems used by the Navy are of the 
same generic type with relatively minor customization except 
in hardware from one installation to the next. These 
designs have been based on analog models intended to provide 
reasonable dynamic responses for pressure, water level, and 
steam flow which are the only state parameters measured by 
conventional controllers. Although these control systems 
have in general proven satisfactory, there is some indica- 
tion [6] that multivariable control systems can effect 
significant improvement in both response time and demands 
on auxiliary equipment. 
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A typical boiler of this type, the Combustion Engineering 
Type V-2M is shown in Figure 2. The main components of the 
boiler are the furnace which is surrounded on all sides 
except the floor by screen tubes, the superheater, main 
riser banks, economizer, water drum, steam drum, desuper- 
heater, downcomers, and screen headers. Navy boilers 
do not include air heaters and the only air preheating 
which occurs is accomplished by passing the incoming air 
over the exhaust stack liner. This procedure generally 
produces an air temperature of about 100-130 °F for an 
ambient of 70°F. In contrast to the boilers used as the 
basis for previous models, the desuperheater is located 
in the water drum rather than in the steam drum. All of 
the steam generated passes through the superheater prior 
to being split out as main superheated steam and auxiliary 
desuperheated steam. 

All of these boilers are of the natural circulation 
type. The "throttle valve" was considered to be in the 
superheated main steam line. 

B. MODEL FLEXIBILITY 

The model was designed to be as flexible as possible 
to permit its relatively simple adaptation to a particular 
boiler. Since the boiler technical manual or equivalent, 

t 

the steam and gas tables, and standard engineering hand- 
books may be presumed to be available, the model was designed 
to utilize data only from these sources for implementation. 
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To convert this input data into the form required by the 
model, an interface program was developed to adjust' 
parameters and determine initial conditions at the operating 
point for which the input data were extracted. This 
approach was taken to avoid the necessity of redesigning 
the model for each operating range or specific boiler 
design being considered. 

C. ASSUMPTIONS 

To permit tractable computation, a number of simplifying 
assumptions were made in the formulation of the model; 

1. Feed water entering the drum mixes thoroughly with 
the drum liquid. 

2. No steam generation occurs in the steam drum except 
through mass transfer at the water-vapor interface. 

3. There is a uniform distribution of heat flux in the 
generating tubes . 

4. The mixture of steam and water in the generating 
tubes is homogeneous. 

5. Quality varies linearly in the generating tubes. 

6. Heat losses to the ambient except through the 
exhaust gases are neglected. 

7. All downcomers and risers in a particular circuit 
have the same length as the average length, and the 
shapes and diameters are equal. Also the number of 
tubes is the same as in the actual boiler. 

8. The slip ratio is 1.0 (no slip). 
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Since the feed enters through a long, horizontal feed 
pipe, and the drum liquid is in a state of turbulence, 
the first assumption is well justified. Although the 
second assumption is not strictly correct during transients 
when there is a significant change in drum pressure, and 
hence saturation temperature, the effect of steam generated 
at the drum wall-liquid interface appears to be negligible 
based on an analysis of this factor in connection with 
this study. The assumption of uniform heat flux is 
recommended by Reference [7] as a reasonable alternative 
to a stepwise integration over the tube height. The 
assumption of a linear change in quality is a consequence 
of the assumption of a linear heat load. For a well 
designed boiler, heat losses through the casing are small 
compared to the total energy generated. The effects of 
any such losses will be reflected in small errors in the 
computed heat capacities of the gas and in the heat transfer 
coefficients. The assumptions regarding tube geometry are 
made primarily for computational simplification. The 
last assumption regarding slip is not an essential one, 
and higher slip ratios could easily be considered if 
necessary to the analysis of a particular boiler design. 

In addition to these explicitly stated assumptions, those 
inherent in the correlations employed in the model formula- 
tion are also considered to have been made. 
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III. MODEL ANALYSIS 



A. GENERAL 

The model was envisioned as an assemblage of control 
volumes for which the rates of energy transfer are 
described by appropriate equations to obtain the rate of 
change of energy storage. The resulting differential 
equations were integrated to obtain the "state" of the 
system. In this way, the states can be chosen to have 
specific physical interpretations. The system is shown 
diagrammatically in Figure 3. 

The following control volumes were considered in the 
model formulation: 



1 . • 


Riser 


tube walls 


2. 


Superheater tube walls 


3. 


Economizer mass 


4. 


Desuperheater tube walls 


5. 


Riser 


fluid (kinetic energy) 


6. 


Riser 


fluid (thermal energy) 


7. 


Steam 


drum water volume 


8. 


Steam 


drum steam volume 


9. 


Water 


drum 


10. 


Furnace 



B. FLUE GAS AND COMBUSTION SYSTEM 

The flue gas and combustion system is shown schematically 
in Figure 4. Heat is carried into the furnace through the 
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sensible heat of the air and fuel, and the heat released 
in combustion. Since the fuel-air ratio is relatively 
constant, and the fuel and air have relatively constant 
inlet temperatures, the sensible heats were lumped with 
the heating value of the fuel. The mass flow rate into 
the furnace is given by 



and the total heat input is 

= Gfc * Ho . 

The heat transfer to the furnace screen is primarily through 
radiation , i . e . , 



The heat balance for the mass of flue gas in the furnace is 



where the last term represents the heat carried off by the 
flue gas. Heat transfer to the remaining tube banks is 
primarily through convection. The heat given up by the 
flue gas is given by the general relation 



GrR. * Cr B +■ G-l 




7 Oii * 4to 



JOC 





^ 3<».s m fA Cp AT 
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and the heat transfer to the tubes may be estimated with 
the Grimson Correlation [8] 






Since the physical properties in the above equation do 
not change significantly over the expected range of 
temperature variations, the correlation is used in the 
form 

‘j-W = K G* CT^-Tw,) 

Since the tube geometry is not known in advance, an n = 0.6 
was used as a reasonable average for most geometries found 
in practice. The equations for the tube banks were: 

1 . Superheater 



ER.J, - Gfc. Caio. (Tr.z “Ta-i) 
EtU, * Go.’^ (Ta&.t -Tiav) 

Tq.o. 2* (Tft-J. * Taj) / X..0 
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2 . Main Bank 



» Gp.Cfc 2 , (»n .3 ” >*m) 

^*•3 * K&j Gc.’^ *“ Trv\a) 

' >7.6.3 ■*• CTl.3 "*■ T&m) j X-0 



3. Economizer 



^K.3 (Tc.m -Tg.s) 
^2-H - Gr!** (Ta<xH — Tw^) 

T«i^ = (Taq + Tas') / X.o 



The arithmetic average temperature for the gas was used in 
the driving potential AT rather than the log mean 
temperature differences (LMTD) , because a comparison of 
the results obtained for the two methods indicated the 
differences were slight for the temperatures involved and 
the simpler average was more computationally efficient. 

C. WATER-STEAM CIRCULATION SYSTEM 

As is shown in Figure 5, the basic circulation system 
consists of the steam drum, water drum or header, and the 
riser and downcomer tubes. The driving force for circulation 
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is the difference in density between the riser and down- 
comer loops. The pressure loss around the circuit balances 
this driving force. For the downcomers, the force balance is 



"Pwi - 



Rw3U*r-(Pj,iirvZ, 







G ujt 3 " 



for the drum downcomers , and 



G ws 2 






2 Ajs 



for the screen downcomer. 

For the risers, the equations are: 



Pwl - Pol = 



+ *“J i (4f 



Pwx - Pot -(F 6 37 ■+• l) 



G?Za: 



+ l c - 



Gu>& 



-h 7 4 



A‘A‘ ’SA s ‘£„ 

+ 3 H„. + ft* £ (£) 



for the screen and main bank, respectively. The friction 
factors are calculated from [9] 




1 

C/.7V - Z toy ) 
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The friction factors for the risers are multiplied by 
Becker's factor [10] to account for two-phase effects 



Be » 1 



/&07g . sg 




The relationship between riser and downcomer mass flow 
rates is given by the continuity equation 

d _ 6w " (>»> 

77 ' av ' v 



The average density in the risers is found from 





From the assumption that the quality varies linearly, the 
density may be obtained as a function of i from 




1 



4 L 



Regrouping the terms we get 



fti) 



1 
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This expression can be integrated in closed form to yield 



- Li 




lx 



- rrHievi, 



Combining terms and simplifying finally yields 



A 






= u 



Xc 



Vj 



+ 1 



From this result, p is obtained as 

' av 



« 



— ks — U 









+• f w Huft 



An energy balance on the riser tube banks yields the 
quality 

rt { <• -T- •++u')} v » 

— WwL Gv (Hwi V EvVr* — Grui (7^u/i **M*v ) 
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where VB is the boiling volume, and the last term on the 
right accounts for heat added in the non-boiling height 
to raise the enthalpy to saturated conditions. 

The non-boiling height is obtained from 

Li - L — 

where q g is the sensible heat added to the liquid in the 

non-boiling height, and q fc is the total heat input. The 

ratio q /q, is given by 
s t 

Hwi. ~ Hwr 
( Hwt +■ X.e “ Huiv* 

The boiling volume is related to the non-boiling height 
by 



Vb 



l» *** w 

u 



With substitution of the appropriate parameters, the above 
equations were applied to both circulation loops. 

The mass balance of the water in the steam drum is 

“ 7 77 = Gro 0 “ * 0 ^ v +• Gb3 t G WX, — Gu>r ~Cru> $ 
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and the energy balance is given by 



(Hvo 1 fowl) * Go (l-xJ Hwi + Gs C'”**) Hwa. 

V Gd 3 W +■ Gwi — Gt»r Unix — Gu»t |4ivX 



Since all the incoming feed must be passed through the 
riser tubes to be heated, an energy balance on the 
downcomers gives 



G to 14 Gwv* t Qru > 5 ■" Gvor 14wi 

— (_ G-ui\r +■ Gu>s) Huix 

Gwt “ \4wi ) 4- (6u»r + ) l4vot 

— ^ QvJr t Hu»t 



Solving for HW2 , the downcomer enthalpy, we obtain 



IW = Sarkis t 14 "” 1 - H “*) ^ h *- 1 
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D. WATER-STEAM HEAT TRANSFER 



1. Risers 

The heat transferred to the riser tubers is used to 
bring the incoming flow to saturation conditions in the 
non-boiling height, and then generate steam in the remain- 
ing portion of the tube. To permit a tractable equation 
set, however, the tube walls were assumed to be at a uniform 
temperature, and all the heat was assumed to have been 
transferred under fully developed flow boiling conditions. 
Since the non-boiling height for the main bank was found to 
be small, the assumption was reasonable for that circuit. 

The non-boiling heights were higher in the screen circuit; 
however, the error in temperature introduced was less than 
three percent which was deemed an acceptable level. For 
fully developed flow boiling, Tong [11] recommends 




A u l3 
ft TsoA 

I.T8H * 



A TV 3 , 



Therefore the heat transfer to the risers was given by 

_ . 4/3 . . 3 

Ew 3 = t Q — / 01 / 

Poi^ Ct« 3 - Tot) 3 



The tube metal temperatures were obtained from an 
energy balance/ 
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Eftl. — £u>3 






<TTmi 

it 



for the screen, and for the main bank 

IAwsCm-^- ERi- e “M 



2 . Desuperheater 

The Dittus-Boelter correlation [12] was used to 
determine the heat transfer from the steam to the tubes, i.e. 

-8 

£d>. 53 ^oxCroo. CtlATb) 

The log mean temperature difference was 

\_MTD = ~ Zkil 

In f-Cfo- 

l_(Tb3 -Tms)J 

The energy given up by the steam was obtained from 

Eox * CdxG-oi. (Toi -To 3 ) 

The heat transfer to the water in the drum was 

Euji> — (.Tms* - ~Tuji ^ 
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The metal temperature was obtained from the energy balance 



lAws Cm 






E.&a - £*>>. 



3 . Economizer 

The Dittus-Boelter correlation gives the heat 
transfer from the tubes to the water as 

£vOi — G-voj. LmtL 



LTV\TQ - 



T*g ^ S 

||J- 

- |w 1 J 



The heat absorbed by the water was 
Evul - G-w^x. Cw ~ TkiJ ^ 

4 . Superheater 

The heat transfer to the steam was given by 

EDI * Kfe-L Gm .' 8 CTfti-TbO 

and the heat transferred from the metal in the tubes 
was 
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r 8 

b-Ot L.MT 0 



liy\tc> 



*Ti bx — Tb*. 

, r i 



The metal temperature was found from 



Eej.- £t.i 

E. PRESSURE DROP CALCULATIONS 

For turbulent tube flow, the pressure drop is given 



by 



A"? «• V 



Dm 



Now 



U m - 



Gr 



) 



SO 



AP 




Lumping the constants we get 
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A? - C 



gL 

Therefore the pressure drops in the superheater and 
desuperheater were given by 

Tit-Wx- c n 4^ 



and 



"Pox - A 3 . c, 

}baX 

The average densities were found from the ideal gas law 
under the assumption that the temperature variation can- 
be neglected, i.e., density is a function of pressure only. 
This assumption is motivated solely by the need for compu- 
tational efficiency. It is reasonable for the desuperheater 
since the temperature excursions there are small. Although 
errors are larger for the superheater density values, they 
are still within accuracy overall limits for the model. 

The flow through the throttle value was given by 

Gr&4 - C 3 "Pfcx 
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F. EQUATIONS OF STATE 

The following equations, obtained from Reference 113] , 
were used to determine state properties. Within the range 
of pressures for which the model is valid (300-1500 psi) , 
the results were within a few percent of tabular values in 
the worst case. 



l>sa t * C \a VW - *MUlo&V .ZIsHSl- \ 

= EXP [ .2ZIS1 W\?S-.V V \ 

^Sf 53 (*3-B - o on0i"Ts^ ■+■ l-\3XE" S' 's«-t 
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IV. DETERMINATION OF CONSTANTS AND INITIAL CONDITIONS 



A. GENERAL 

For the model to be useful, it should be possible to 
obtain the necessary constants and initial conditions from 
the data supplied in the boiler Technical Manual augmented 
with reasonable engineering estimates and tabular data from 
standard handbooks. 

The primary items of interest are heat transfer coeffi- 
cients and initial conditions such as steam flow rates , inlet 
and outlet temperatures, circulation rates, etc. While 
some of these were directly obtainable from the boiler 
Technical Manual, others such as tube wall temperatures 
and riser exit quality were not. Therefore, a program was 
developed to permit the user to enter known data elements 
from the Technical Manual and other sources, and to obtain 
as an output the necessary values for the constants and 
initial conditions to ensure static balance at the speci- 
fied operating point. 

The basic premise under which the program operates 
is that the model equations with all derivatives set to 
zero form a determinant set with the addition of several 
relations which apply under equilibrium conditions. These 
relations are: 

a. Riser and downcomer mass flow rates in 
each circuit are equal. 
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b. Total steam flow from the boiler equals 
the incoming feed flow rate. 

c. Total heat transfer to the incoming feed 
equals the heat transfer to the screen and the 
main bank, and the heat removed from the steam 
passing through the desuperheater. 

d. The net mass transfer at the steam liquid 
interface in the drum is zero. 

B. PROGRAM EQUATIONS 

The total heat input to the boiler was given by 
^6-leAsfc /sec.-f t 3 ) * Fv>r.»»*c 6. Volume C'fi' 3 ) 

If the sensible heats of the oil and air are lumped with 
the heating value of the fuel, the heating value is 

Bo - Eft /G-6 

The heat absorbed by the furnace screen at steady state 
was calculated from 

F.P.1 - Fufl.HAC.£ BfcAT / sot -£t*) 

X BtCcV & Vv\«j 
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Since the flue gas temperature and the screen tube wall 
temperatures were provided, the heat transfer coefficient 
for radiant heat transfer to the screen could be determined 



from 



fa* - 



Tiru 



/OO 



This result holds for the conditions reported in the 
manufacturers technical manual. Degredation due to age or 
fouling must be accounted for through the use of appropriate 
correlations . 

With the enthalpy of the feed water at the economizer 
outlet and the enthalpy of the steam leaving the dry pipe, 
and the change in enthalpy of the steam flowing through the 
desuperheater, the heat transfer to the main bank was 
calculated from 



The metal temperature of the main bank risers were obtained 
by calculating the constant 



E&3 » G-64. (Mt>i-Ww4 S ) -£ai. "6 m.(H SH owT - H w#wT ) 





nfex x 10 *■ 



and solving the heat transfer equation 
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Tiftj - 



£ Eb .5 { ^ ^ 



'/3 

+*r di 



The specific heat of the flue gas passing the main bank 
was 



CRj, ” / G-h- C~Tn.3 - I b-v/) 



and the unknown constant in the gas side heat transfer 
relation was determined from 

Ke.3 = / Go} -tv*) 

In a similar manner, the specific heat of the furnace 
flue gas is 

Cat - C e a- fcai.') / <Sft.Tai 



and the water side heat transfer coefficient is 



K«J = E<u / Vot 13 (T„t -Tot ) 3 . 

The heat transferred in the desuperheater was: 

E 0r = G&t ( V\s« oar - Hw wT ) 
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With this value and the temperature difference the specific 
heat of the steam was computed from 



^ * Eo%. / Got (Tbt -Tfcs) 



The heat transfer coefficient was calculated from the 
Dittus-Boelter correlation 111] 



Kg'l - 




The log mean temperature difference was 

LWTb ~ Gbl. 



and the metal temperature was determined by 



“Tiwy * Cic>t"Q»od')/C»-c3^ 

where : 

§= <=xp [ 

With the metal temperature, the water side heat transfer 
constant was obtained from 

Ku> 2. * ^Oi / CTi v\s — Tuii 
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Analysis of the Economizer is more complicated since it is 
a multipass heat exchanger with finned tubes. Additionally, 
since we are concerned with energy storage within the metal 
of the economizer, the internal and external heat transfer 
must be considered separately. 

The heat transferred was given by the change in enthalpy 
of the feed water, i.e.. 




The specific heat of the water was 




As in the case of the desuperheater, the inside heat 



transfer coefficient was calculated from the Dittus- 



Boelter correlation as 




From this the LMTD was obtained as 



LttTb* J £u.i 
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Finally the metal temperature was determined from 



- (3Twj ) / ( I -ft) 
P * { ‘■"to /c-r^ --Wi) \ 



The outside heat transfer constant was 

« Etoi / 6(2.’ ** CTftAi -Tuq') 



where 



Tl,j = — jL t* T ^~ 

The specific heat of the flue gas passing the economizer 
was 



= Ew± / 6(2. c 1 M 

The superheater heat transfer was obtained from the change 
in enthalpy of the steam, i.e., 

<s-&-l(;h s „ w . m 01 .) 

With this value and the known metal and gas temperatures, 
the specific heat of steam, the specific heat of the flue 
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gas passing the superheater, and the inside and outside 
heat transfer constants were determined as 



Cot = Efel/Obt CTot-Toi) 
£oi/ 6 a c Tat -Taj) 

i<0 ' L = L ^ r , b 

— / &l 2. ^ C T“/2A1. — J ryix) 



where: 



LWTO = 



TVx ~ 'T'ryi. 



and 



T«Ai - CT(U 4 Till) / X o 

The superheater and desuperheater pressure drop constants 
were evaluated from the known equilibrium pressures as 

£*1 = C ^ / £ 0 i 1 

^ - (.Pox “”^3 ) Ph«.X /c DX Z 
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Similarly, the constants used in the calculation of the 
average steam densities were obtained from: 



C ** /(fo t 



and 



c, 



10 * ^ / c a t * ^ ) 



Finally the steam value flow constant was computed from 



C. CIRCULATION BALANCE 

Since they are not normally provided by the manufacturer, 
values for the friction factors, and bend and entrance 
loss coefficients must be assumed. The values do not appear 
to be particularly critical, and tables found in engineering 
handbooks may be used to determine reasonable values for 
these parameters. 

An iterative procedure must be used to adjust the various 
flow values involved to achieve a steady state balance. 

Since the riser exit quality for boilers of this type 
generally lies in the range of .02 to .04, a starting trial 
value for the circulation flow rate may be obtained from: 
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E.* = fry. 

where X E is the exit quality, E^ is the heat transfer to 
the riser, G r is the flow rate, and HFG is the heat of 
vaporization. With this trial value, the following 
balancing procedure is used: 

1. Set HW2 (Downcomer Enthalpy) = HW1 

2. Calculate HWG = HW2 + ED2/G0 

3. Calculate the riser exit quality: 

^ = 4- VAu>i — ^ Gr a 

4. Determine the non-boiling height: 

Wn&ft. « Uxfc Ct^iot - / C + Xo — Hww) 

5. Compute the average riser density: 

6. Obtain the outlet density: 

* l o / CV4 
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7. 



Compute the friction factor: 



Fa - * { I.m aiu) I l 1 * W'O’as* 

8. Next the riser loss factor is computed: 

*&/?*«. * Zh/^ 

9. Finally we may calculate G0': 

7 Vi. 

Gcf - ^ (.£*» ^ / C Ctoo +■ ^0 i 

where Cl 00 is the downcomer loss factor. 

10. Compare G0 and G0 ' if equal go on, otherwise 
update G0: 

Go - - £<?'“&* + 
and return to step 2. 

11. Assume a trial value for Gs (first pass only) . 

12. Compute Xs : 

)(<* * Co-oi - Go Xo ) /6*s 
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13-19 Same as for the main bank risers. Return to 
12 if reiteration necessary. 

20. Obtain HW2 ' from: 



Hwj. — — Mu-*) +* 



21. If HW2 ' = HW2 go on, otherwise update HW2 ; 



Mwi - + 

z 



and return to step 2. 

Upon exit from the iteration loop, equilibrium values 
of the flow parameters will have been established. With 
these and the previously calculated heat transfer parameters, 
the program generates the list of constants and integration 
initial conditions required by the dynamic model. The 
program listing and a sample output is given in Appendix B. 
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V. RESULTS AND CONCLUSIONS 



A. OPEN LOOP RESPONSES 

The output curves for 10% increases and decreases in 
throttle valve opening and fuel flow are given in Figures 6- 
33. Ramp inputs of approximately 1.5% per second were used 
to simulate flexibility test conditions. Only one input 
parameter was varied for each run with all others held 
constant. The results were similar to those presented in 
[3] and [4] insofar as general trends were concerned. Since 
the boilers involved were of widely different geometries 
and sizes, a quantitative comparison was not possible. 

Since CSMP uses only single precision arithmetic packages 
and integration algorithms, the high frequency noise in the 
outputs is attributable to the roundoff error inherent in 
the wide variations in magnitudes of numbers used in the 
model, e.g., quality (order of 0.01) and riser mass flow 
rate (order of 1000) . In common with the results of [3] 
and [4] , the "swell effect" in response to an increase in 
throttle opening is negligible. It should be noted that the 
swell exhibited in the data of Whalley [3] is more apparent 
than real, inasmuch as the scales in his graphs are very 
magnified and the rise is actually only 0.01 inches. The 
sudden increase in riser mass flow rate shown in [3] is 
also not present in the current model. 
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Since "shrink" and "swell" effects are of some concern 



to controls designers, an attempt was made to model for 
such effects. Four possible causes for the phenomenon 
were considered. 

1. Homogeneous Nucleation in the Steam Drum Liquid 
This effect is caused by a difference between the 

saturation temperature corresponding to the actual drum 
pressure, and the actual liquid temperature. The super- 
heat required for such nucleation is quite high however and 
nucleation will occur at the heating surfaces at much lower 
superheats. Therefore this was considered to be an unlikely 
cause for the effect. 

2 . Nucleate Boiling at the Drum Wall- liquid Interface 
The driving force for this phenomenon is also the 

temperature difference between the drum surface and the 
saturation temperature. Since the temperature difference 
will be slight except under large drum pressure excursions, 
the impact of this factor should be slight. The Rohsenow 
correlation [5] may be solved in terms of a rate of vapor 
generation •• to yield 




ci l 3 / H& Ba 

c» fv / 1 ] v ^ tir 



MiH 



It was assumed that the vapor formed would be uniformly 
distributed bubbles which would rise to the surface with a 
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velocity given by [14] 



°\> * iV7§~b~ 

which is a conservative estimate since interaction effects 
tend to increase the rise velocity. The bubble diameter 
is given by [14] 



X) = 24* ? 



This gives as .25 ft/sec. The average liquid depth is 
15 inches, so the transit time is = 5 seconds. The effect 
was modeled crudely by assuming the entrained bubble 
mass was 5.0(q/hfg) lbm. The impact of this model on shrink 
and swell values was negligible, since no large pressure 
excursions were developed. 

3 . Expansion and Contraction of Varpor Bubbles Entrained 
in the Riser Fluid " 

An attempt was made to model this effect by including 
the riser tubes in the control volume used as the basis for 
liquid level determination. The equation used was 

Vol = * (■&« ~ 1 ) Ae 

where Ralx and Ra2x are the average riser densities computed 
using the specific volume for the saturated vapor corresponding 
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to drum pressure, while Ral and Ra2 were computed from 
the liquid temperature saturation conditions . The 
effects were insignificant. 

4 . Transient Deviation of Riser and Downcomer Mass 
Flow Rates 

The small swell effect obtained by Whalley [3] 
can be attributed to the temporary increase of riser mass 
flow rate over downcomer mass flow rate. The period over 
which this imbalance could be sustained is short and would 
fail to account for the longer duration of "shrink" and 
"swell" effects observed in practice. It may, however 
trigger one of the two phase flow instabilities common to 
natural circulation loops e.g., flow oscillations [10,13]. 
The analysis of these effects was not exhaustive, and the 
modeling attempts were crude. In the case of the first two, 
however, the contribution to swell is probably not signifi- 
cant. Modeling difficulties with the latter two effects 
precluded any specific conclusions. 

B. CONCLUSIONS 

As with all simulations , the information available from 
the model is extensive. Notably, in this study, main bank 
and water-screen circulation values can be separately 
evaluated. Hence, the contribution of each of these ele- 
ments to overall boiler operational response and efficiency 
can be determined. 

The model simulates the responses of a Naval boiler 
quite well. Unfortunately, open loop data is unavailable 
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for the actual boiler and insufficient information was 
available to develop a model of the existing LHA-1 control 
system as a basis for comparison studies. As such, however, 
the model forms a basis upon which future studies involving 
modern control system designs can be established. Moreover, 
while considerable refinement of the model is possible, 
the advantages of greater accuracy must always be weighed 
against increased computation costs. In this regard, the 
implementation of more efficient integration algorithms 
would be particularly beneficial. 



55 



< feed in 



G01 



TW3 


ECONOMIZER 


HW3 






HW4 




TW4 



DRUM 

downcomer/ 



/ 






G03 

iU 



HOI 

VD1 

PO 

TO 






STEAM 

v/ORUM 



V 



- Hwl 



7 



TDl 



POl 



MAIN 

BANK 



Ral , 



lO 



WATER 
" 0 R UM 



Hw6 



VG04 
PD2 
T 02 



G02 
P03 
T 03 



THROTTLE 

VUVE 

-tJj— 



v Gs 



HEADER 

DOWNCOMER 



Xo 

Ro 




-C>MAIN STEAM 



— P> A UX STEAM 



FIGURE 1. SIMPLIFIED BOILER SCHEMATIC 
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FIGURE 2. COMBUSTION ENGINEERING TYPE V-2M BOILER 
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FIGURE 3. CONTROL VOLUME ARRANGEMENT 
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FLUE GAS AND COMBUSTION SCHEMATIC 
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FIGURE 5. BASIC CIRCULATION SYSTEM 
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FIGURE 6. WATER LEVEL RESPONSE TO A 10% INCREASE IN THROTTLE SETTING 
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FIGURE 10. SCREEN CIRCULATION RESPONSE TO A 10% INCREASE IN THROTTLE 
SETTING 
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FIGURE 12. AVERAGE RISER DENSITY RESPONSE TO A 10% INCREASE IN THROTTLE 
SETTING 
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FIGURE 13. WATER LEVEL RESPONSE TO A 10% DECREASE IN THROTTLE SETTING 




69 



r 0.00 30.00 60.00 90.00 120.00 150.00 160.00 2'l0.00 290.00 2*70. 00 300.00 

TIME (sec) 

FIGURE 14. DRUM PRESSURE RESPONSE TO A 1Q% DECREASE IN THROTTLE SETTING 
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FIGURE 17. SCREEN CIRCULATION RESPONSE TO A 10% DECREASE IN THROTTLE 
SETTING 
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FIGURE 18. RISER EXIT QUALITY RESPONSE TO A 10% DECREASE IN THROTTLE 
SETTING 
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FIGURE 20. WATER LEVEL RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 21. DRUM PRESSURE RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 22. STEAM FLOW RATE RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 23. MAIN BANK CIRCULATION RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 24. SCREEN CIRCULATION RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 25. RISER EXIT QUALITY RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 26. AVERAGE RISER DENSITY RESPONSE TO A 10% INCREASE IN FUEL FLOW 
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FIGURE 27. WATER LEVEL RESPONSE TO A 10% DECREASE IN FUEL FLOW 
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T I ME (sec) 

FIGURE 28. DRUM PRESSURE RESPONSE TO A 10% DECREASE IN FUEL FLOW 
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TIME (sec') 

FIGURE 29. STEAM FLOW RATE RESPONSE TO A 10% DECREASE IN FUEL FLOW 
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FIGURE 30. MAIN BANK CIRCULATION RESPONSE TO A 10% DECREASE IN FUEL FLOW 
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TIME (.sec) 

FIGURE 32. RISER EXIT QUALITY RESPONSE TO A 10% DECREASE IN FUEL FLOW 
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FIGURE 33. AVERAGE RISER DENSITY RESPONSE TO A 10% DECREASE IN FUEL FLOW 



APPENDIX A 



ThI3 PROGRAM CALCULATES THE CONSTANTS AND INITIAL CjNCITICNS 
FCR A US NAVY D-TYPE BCILEF GIVEN THE INPUT OA T A RETAINFC FRO” 
ThE ECILEP TECHNICAL MANUAL * STEAM TA6L : S, AMC STANOARC SNGINE.-- 
ING HANDBOOKS. THE INPUT CAT A ARE ENT ERGO IN NAVE LIST FORMAT. 
NAMELIST DESIGNATIONS ANO INPUT N-HCNIC' FOLLOW 



INCCM: 

OPPNT 

TCTSTM 

SPHSTM 

CShSTM 

PC« 

ShCT 

ShCP 

DShCT 

OSbCF 

= CCN IT 

z C CNOT 

AISTMP 

1 I LTMP 

A I3FLC 

11LFLQ 

XCSAIP 

DR4F7 

TG* SSC 

TC^SSb 

TG A SMB 

TGASSC 

TSCRN 

TSFbTR 

bRFVQL 

FhAS 


OPERATING POINT (PERCENT ) * 

TOTAL STEAM F LOW ( L B/HR ) 

SUPERHEATED S TEAM FLOW (LB/FR) 

Of SUPERHEATED STEAM FLOW (LE/HP) 

DRUM PFESSUR C (PSIG) 

SUPERHEATER CUTLET TSMPERATlR c (COO-F) 

SUPERHEATER OUTLET PRESSLR C (PSIG) 

CESUF-PHEATER CUTLET TEMPERA-UP £ (CEG-F) 
DESUFFFHEATEP CUTLET PRESSURE (PSIG) 

ECONCMI Z C P. FEED INLET TEMPERATURE (OcG- c l 
ECONOMIZER FE C C OUTLET TEMPERATURE (C C C--F) 

air inlet temperature to registers (cec-f) 

OIL TEMPERATURE aT BURNER INLET ( OEG-F ) 

AIRFLOW R AT c (Le/hR) 

CIL FLOW R A 7 C (LB/HR) 

EXCESS A Ip (PERCENT) * 

DRAFT LOSS (IN-H20) 

FLOP GAS T EMP C RATURE LEAVING SCREEN ( C r G- F ) 

RLU: GAS T=mp=catj*s LEAVING SUPERHEATER (CEG-F) 
FLUE GAS TEMPERATURE LEAVING MAIN BANK (OEG-F ) 

FLUE GAS TEMPERATURE LEAVING ECONOMIZER (OiEG-F ) 
SCREEN TUBE WALL TEMPERATURE ( DEG- F ) 

SUPE FHEA T EP. TUBE WALL TEMPERATURE (CEG-F) 

FEAT RELEASE (KSTU/HR/CU FT C UPNACE VOLUME ) 

FURNACE HE A T ABSORBT ICN (KETU/HP/SG F T RADIANT HR AT 
ABSORBING APEA) 



NOTE: "PERCENT" VALLES ARE TO BE LOGGED AS DECIMALS. 



iSCRcEN 

DTIBSC 
L A V SC 
NTLBSC 
Rb tSSC 
MASSSC 


AVERAGE. INSIDE DIAMETER OF SCREEN TUBES (IN) 

AVERAGE LENGTH OF SCREEN TLE5 (FT) 

NUMBFP OF SCREEN TUBES 

RADIANT HEAT ABSORBING AREA OF FURNACE SCREEN (SC FT) 
TOTAL WEIGHT OR SCREEN TUBES (LB) 


$ S PHTR 

AS? ASH 
N* A S S SH 

$ N E A NK 


TOTAL SUPERHEATER HEAT TRANSFER ARFA (SC FT) 
TOTAL WEIGHT OF SUPERHEATER TUBES 


OTU0M8 

LA 

NTLBMB 
NASSMB 
A R c AM B 


AVERAGE INSIDE DIAMETER OF MAIN BANK TUBES (IN) 
AVERAGE LENGTH OF MAIN BANK TUBE (IN) 

NUMBER OG MAIN BANK TUBES 

TOTAL WEIGHT CG MAIN BANK TUBES (LB) 

TOTAL HEAT TRANSFER AREA (SC FT) 


$ 5 C C N 

0TL3RC 

NPASSE 

NTLBSC 

LTLBGC 

MASSF.C 


INSIDE CIAMFTER OF ECONOMIZER TUBES(IM) 
NUMBER C e TijPE PASSES 
NUMBER OF TUBES PER PASS 
AVERAGE LENGTH OF ECJNJMI2ER TUBE (F^) 
TCTAL WEIGHT OF ECCNCMIZER (LB) 



$ C £ SP F 



nnnnnnnnnnnnnnnonnnniinnnnnonnonnnnnnnnnnnnnonnrtnnnononnnnnnonnonnnonn 



CTL2DS INSIDE D I A NOTE P CF OcSUPEPhS-ATER TueE (IN) 

NTIBCS NUMBER OF CESU FERH :A TER TUBES PER PASS 
NFASSC NUMBER OF CFSU FERHzAT r.R TIE' “ASSES 
LTL3DS LENGTH OF CESUPERH 2 a T ER TIES (FT) 

AREACS TOTAL HEAT TRANSFER AREA CF DSSUPE PHc AT“R (SC FT) 
MASSDS TCTAL WEIGHT CF DESUPERHEATER ASSEMBLY (LB) 



SCFMCCP 

CTLECC AVERAGE 0 I AMETER OF DRUM CCWNCOM-R TL8ES (IN) 
LAVDD AVERAGE LENGTH OF DRUM DOWNCOMER T UE'' (FT) 
NTLBCD NUMBER OF CRUM DOWNCOMER TLEES 



SHCRCCR 

CTLEHC AVERAGE INSIDE DIAMETER CF HEADER CCWNCO 1.ER (IN) 
L A V HD AVERAGE LENGTH CF HEADER CCWNCOMF.R (FT) 

NTL3HD NUMBER OF HEADER OCWNCOMERS 



SECILER 



CSTNDN DIAMETER OF STEAM DRUM (IN) 

LSTMDM LENGTH O c STEAK DRUM (FT) 

CWTFCM Cl AM FT FR OF VIATOR DRUM (IN) 

LWTPDM LENGTH OF WATER DRUM (FT) 

HNCRM HEIGHT OF NORMAL WATER LEVEL ABOVE BENCH MARK (FT) 
HHC<? HEIGHT OF HEADER (SCREEN) AeOVE BENCH MARK (f=T) 

HWTRDM HEIGHT OF WATER DRUM ABOVE BENCH NA°I< (FT) 

FL P VOL FURNACE VOLUME (CU FT) 

NCTE: CHOICE OF B C NC H MARK IS ARBITRARY 

$ THE F M C 



HSHCUT ECONOMIZER GU T LET ENTHALPY (BTIJ/LBM) 

HD SCUT DESUPERHEATER CUTLET ENTHALPY (BTU/LEM) 

HECIN ENTHALPY OF ECONOMIZER FFEC INLET ( B T L 1 / L B M ) 

He C OUT ENTHALPY OF ECONOMIZER FE EC OUTLET (BTL/LBM) 

(1IKH2C THERMAL CONDUCTIVITY CF WATER ( 9 T U/ F R . F" . C E G - F ) 

(DPRH2C PRANCTL NUMBER F JR WATER 

( 1 ) VSCR2CI KINEMATIC VISCOSITY FDR WATER (L3M/FT-SEC) 

(Z)KSTM THERMAL CONDUCTIVITY FOR STSAM ( BTU/HP .F?. C2C-- F) 

(2IPRSTM PRANCTL NUMBER FOR STEAM 

( 2 ) VSC STM KINEMATIC VISCCSITY FOR STEAM (LBn/FT-SEC) 

RSHOUT SUPiPHEATER CUTLET DENSITY (LBM/CU FT) 

PCSCUT CESUPrRHEATFR CUTLET DENSITY (LBM/CL FT) 

RFLLE OENSITY OF FLUE GAS AT TGASSC (LBM/CL FT) 

NOTES: (1) FVALUATEO AT AVERAGE ECONOMIZER WATiiR TSNFcRAT|_P" 

(2) fcVALU AT FO AT AVERAGE CE SUPERHEAT SR STEAM TEMPERATURE 



KSCSC ROUGHNESS RATION (SAND ECU I V ALEN’’ ) FCR SCREEN TUB-'S 

KSCMB ROUGHNESS PATIO (SANO EQUIVALENT) CF MAIN BANK TUE' 

KSCCD ROUGHNESS PATIO (SANC EQUIVALENT) F C F DRUM CCWNCCM- 

KSCHD ROUGHNESS RATIO (SAND EQUIVALENT) FC F HEADER CCWNCC 

ENTSC SCREEN ’•(JPE ENTRANCE LOSS 

BEN CSC SCREEN TUe? BEND LOSS 

HNTMB MAIN BANK ENTRANCE LOSS 

BENCMB MAIN BANK BEND LOSS FACTCR 

ENTCC CRUM DOWNCCMER ENTRANCE LCSS F AC T OR 

BENCCD DRUM DCwNCCMER BEND LCSS FACtqp. 

ENTHD HEADER DOWNCOMER ENTRANCE LGSS <=AC t OP 
BENOHD HEADER DOWNCCMER dSNC LOSS FACTOR 



IMPLICIT REAL(A-Z) 
INTEGER I 
INTEGER J,K,M,N 
INTEGER NPASSc 
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c 

c 

c 



DIMENSION TH20( 20 ) t CPASS (20) 

DEFINE NAM'S L I ST 5 

NAMELIST/TNCCNO/CPPNT ,TCTSTM, SPHSTM, CShSTM,P0PUM ,SHCT ,SHCP, DSHCT , 
lDShQP,2CONIT,ECCNOT,AlPTMF,OILFLO, A I R F LO , XC S A I R , C P A FT , T G A SS C , 
2TGASSH ,T GAS MB, TGASEC. TSCRN,TS PH TR, FRF V CL t P HA S , C I LT M P 
NAMELI ST/SC RE EN/D?UB$C,L A VSC, NTUBSG , P F ASSG , M AS S S C 
NAMEL ! ST/S D hTR/AP FASH, MAS SSH 

N A M E L I S T / M f I 8 A N K / CT U 3 M B , L A V M B , NT U 3 M B , M A SSM 6 , ARE A M E 
NAMELIST /EC ON /DTL'BEC , N PA S SR , NTU B EG , LTL E EG , M AS S c C 
NAVfL 1ST /C c SP H/ CTUB OS tNTUBDSfNP AS SO f L1L3DS,AF r AC S, MAS SCS 
NAMEL IST/CRM0GP/CTU3CC tLAVDOt NTUBOC 
NAM EL 1ST/FDRDCR/DTUBHD ,LA VHD, NTU6HD 

NAMEL IST/BOILER/CSTMICM,LSTMDM,DWTRDM,LWTRDM,HNCFM ,FHDF ,FWTRDM , 

$F L R VC L 

NAMEL IST/THERMO/HSHOUT,HDSOUT,hSC IN ,FEC3L'T, 
lKF20,PRH2CfVSCH2C,KSTM,PPSTM, VS CSTM , R S FOU T , ROS CLT , RF LL E 
NAMELIST/LOSSES/KSDSC , K SC MB , K SDDO , K S C F C , EN T S C , B l N CS Z , E NT M E , 

$BEN DM S , 5NTDC , BENC0D » r NTHD , 8LNDHC 
NAMELIST/ IN CONI /TRIO, TM13 ,tm 23, TM30 ,TM 40,TM50 , C-CRC ,GSRC ,C >CC ,Q /SC 
NA^E LIST/ INCGN2 / HM10 , MW 10 t VRQ 3 , NW6 0 ,GW2, TW3 , GL , GE ,GC2 , GC10, V : NT 
NAMEL 3 ST/ IN CONE /DR A 10, CR A 20 , GDI E , R A 1 X C ,RA2X3,GWFC ,GwSC 
NA.v=LIST/C3NSTl/KR,KS ,KCR ,KDS,4R,LP,CP ,LR1, AS, LS, CS,LS1,LCP, AD* 

NAM EL IST/C0NST2/0DRfHNBRt CDSf LOSiHNBS ,iQS,MM , N M 2 ,MM3, MM4,MiM5 ,FXR 
NAME LI ST/C DNS T3/FXS,M>F1, VR,VS ,VNDP,VCP ,LQ1R, LC1S 

NAMELIST/C0NSTA/GR1 ,CRlA t CR2t CR3 tCDI ,CC2,CW,CM,KC1,KD2,KF1,KC2,KR2 
NAMELIST/CDNST5/KR4,Z1,Z2,Z3, Z4,Z5,Z6, Z7,Z8,HW2 , C*C i 3 , G Z2 , C8 , C 9 , C 10 
NAMEL IST/CQNST6/G11 , G7 , CS L , KW 1 , Kw 2 , *W E t KW 4, HC 

N AM EL I ST/0UTPLT/ER1 ,E F2 , E R3 ,ER4 , - 0 2 , GO , GS , X 3 , X S , R A1 ,RA2 , ROt FS , 
$RWi,FWl,FW2,FW2,FW4,HW6,HFG 

READ ( 5, ING3N3) 

WRITE (6,INCCN0) 

READ ( 5 , SCREEN) 

WRIT5(6t 601 ) 

WFITE (6, SCREEN) 

R E AD ( 5 t SPHT R ) 

W F IT E ( 6 , 60 1 ) 

WR ITE ( 6 t S PHTR ) 

R t AC ( 5 t MNBANK ) 

WP ITE (6,601 ) 

WRITE (6 t MNBANK) 

R EAC ( 5 t EC ON ) 

W 9 IT E ( 6 , 60 1 ) 

WRITE ( 6 i EGON ) 

REAC( EtDESPH) 

WFITE (6,601 ) 

WRITEC 6.0ESPH) 

Rt AC ( 5 , CRMDCR ) 

W P IT E ( 6 , 601 ) 

W P I T E ( 6 , DRM DC R ) 

READ ( 5 , FDRQCR ) 

WFITE (6 ,601 ) 

WR ITE ( 6 , HDRCCR ) 

READ (5, BOILER) 

WR ITE (6,601 ) 

WR ITE (6 ,BCI LER) 

REAC(5, THERMO) 

WRITE(6,601 ) 

WRITE(6, THERMO) 

R£AC(E, LOSSES) 

WRITE (6,601) 

WP ITE ( 6,LCSSES) 

CALCULATION CF PRELIMINARY CONSTANTS 

P I -3 • 1 4 1592 7 
G = 2 2 . 2 
CM*. II 
D R A 10 = 0 .0 
D F A2 C - 0 *0 
ZER3 = C .C 
VENT=CPPNT 
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KR^KSCRB 

KS = K$ CSC 

KCR = K SCDO 

KCS=KSDHD 

Z 1 = EN TC C 

Z2=BEHCCD 

Z 3 = cN 7N B 

Z4 = BE!^CMB 

Z5 = SN T$ C 

Z 6 -8 E N C SC 

Z 7 = zN 7 H C 

Z 6 =B E H Q HD 

P C 1 = P CR IM + 1 4 . 7 

H k 1=EXP (0.264 52 *AL0G( PCI) <4.467C6) 

HFG=922.15-C.40516*F01*1.717E-04*PCl**2-4.219E-C8*P01*+3 
HD 1= H h 1 + HFG 

T C 1 = F X P ( .22 151 *AL0G( PCI) +4.7712 3) 

GB =G I LFLC/ 3600. C 

GL-A IRFLC/360G.C 

TR13=7GASSC 

TN 10=7SCRN 

TM 20 = 7 SPhTR 

7R2=TR10 

7R 2 = 7 CASSH 

TR4 = 7 GASMB 

7R5 = TC-ASEC 

GC 1=TCTS 7 M/3603.0 

GC1£=CC1*.3001 

GD4 = S FHS7M/2600.0 

TC2=ShCT 

P D2= S HC P+ 14 . 7 

TC 2 = D 5HC7 

PC2 = CSHCP+14 .7 

GC2 = CSHSTH/3601 .0 

TV\2 = EC0NI T 

TU4 = ECC,NCT 

E B = HR F V0l*F UR V CL* 1000. C/ 3600.0 
ER 1=FH AS*RHASSC* 1000. C/ 2600.0 
L S 1 = L £ V SC. 

LS=MIBSC*L AVSC 

C S = CJLB SC TUB SC/12. C 

AS = FI *CTlJBSC**2*NTUBSC/< 4 .0*144.0) 

VS=4S*LAVSC 

KK1=MASSSC 

M2 = ^ASSSF 

L R 1 = L ^ VMB 

LR=N71BM8*LAVM3 

0R=D71EM6*N TUB ME/12 .0 

AR=PI*0TU&MB**2*NTUB^e/(4.0*144.0) 

VR=AR*LAVMB 

^M2 = M$SMB 

MM4=MA5SEC 

HM5=HASSDS 

LCR=LAVDD*NTUBOC 

LD S = L A VHD*NTUBHC 

CCR=CTU8D0*NTUBC0/12.C 

0CS = D7LBHD*MUBHD/12. C 

ACR = NTUBDD* (P I* (CTUBDC /24 .0 )**2 ) 

ACS = N71EHC* <P I* ( CTUBHC/2 4 .0 )**2 ) 

V D R= ( PI*OSTKDM**2/(4.0*144) ) * l S T M 0 N 

V N CR = (PI*CWTRCM**2/ (4.0*144) ) *1 WTROM 
Hr»3 = HEC IN 

H k 4 = H E CCLT 
GR=GB<GL 

Rfc=62.8-0.01781*TD1+1. 132E-05 *T Cl**2- 6 .786E-G8 *TC 1**3 
VF=1/R h 

VFC=E24.0/PC1-. 1 
RD 1 = 1 • C / ( V F + VFG ) 

R C 2 = R $ H CUT 
R C 3= R C $ OUT 
7 V\ 1 =T C 1 
G V> 2 = G C 1 

H P 1 = P F L LE*FUR VOL 
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CZ2=(CSTMDM/12.C)*LSTMDM* (1.0/12.0) 

CALCULATE LOWER HEATING VALUE OF FUEL CIL 
H0=EB/GB 

COMPUTE STOICHIOMETRIC AIP/FUEL RATIC 
C SL= ( 1 .0-XCSA IR )*GL/G8 
RADIATION HEAT TRANSFER COEFFICIENT 

KF1=EF1/ ((( TRIO -*460.0 ) / 100 .0 ) **4- ( ( TM 1 C + 460 . 0 ) / 10 C . C 1 * *4 ) 

CCMPUTS HEAT TRANSFER TO Th= MAIN TUBE BANK 
ER2=G:i*(HDl-HW4 )-GD2* (HSFOUT-HDSOUT )-ERl 

CONPUTE MAIN BANK AVERAGE METAL TEMPERATURE 

KM S (1.C/1. 782EC6)*ARFAME 

TM2T=(ER3/( KW4*FD1**( 4.0/3.0J)J**(l . C / 2 . 0 ) +TD 1 

C COMPUTE MAIN BANK GAS SICE HEAT TRANSFER COEFFICIENT 

K F 2= E P 2/ ( (( TR2+TR4) /2 .C-TM30)*GR**.6 ) 

C CCMPLTE SCREEN EANK WATER SIDE h£ AT TRANSFER COEFFICIENT 

KW2=ERl/( (PCI** <4. 0/3.0 )*(TM10-T01 J**3 J 
C HEAT TRANSFER TC THE SUPERHEATER 

6R2=GC1*(HSHCUT-HD1 ) 

C 

C SUPERHEATER GAS SIDE HEAT TRANSFER COEFFICIENT 

KR2 = ER2/( (( TR2 + TP3)/2.C-TM20)- i ‘GR 1) ‘*.6) 

C HEAT TRANSFER TC THE ECONOMIZER 

ER4=CC1*( HW4-HW2) 

C COMPUTE SPECIFIC HEATS OF FLUE GAS 

C FURNACE 

CR1= (EE-ER1 )/ ( (TP10-8C.0 MGR) 

C SUFERHEA TER 

CF1A=£R2/(GP*(TP2-TR3 ) ) 

CC1=EF2/(C-D1*(TC2-TD1 ) ) 

C MAIN BANK 

• CP2=ER3/(GR*(TR2— TR4J 1 
C ECONOMIZER 

CF2 = ER4/(GR*( TR4-TR 5 ) ) 

C SPECIFIC HEAT FOR WATER 

CN=EF4/(GD1*(TW4-TW3) ) 

C COMPUTE AVERAGE HEAT TRANSFER COEFFICIENT 

CEC1=CTUE2C/12.C 

AREAEC=NTUBEC*LTUEEC*PI*C5C1*NPASSE 

KWl=AREAECMKH2C/DeCl*n.8)*.023*(4.C/(PI*VSCH2C*NTUBEC ) )** . 8 * 
$PPH2C** .4 

LMICEC = cR4/(KW1 3 *GD1**.E) 

XXXX1 1=EXP( (TW4-TW3) /LMTCEC) 

T M40 ® ( T W3-T W4*XXXX11) / ( 1 . O-XXXX 1 1 ) 

KP4=£R4/(GR**.6*( (TR4+TR5)/2. 0-TV40) ) 

C SPECIFIC HEAT OF STEAM AT AVERAGE CO-SUPERHEATER TEMPERATURE 

ED2=(GC2* (HSHOUT-HDSCLT) ) 

CC2=EC2/(GD2*(SHOT-DSHfT) ) 

C DESUPERHEAT ER STEAM SICE HEAT TRANSFER COEFFICIENT 

DDSH1-DTUBD S/12.0 

H X F R C S = (KSTM./DDSHl**1.8) * .023* ( 4. 0*GC 2 / ( P I*VSC STM*NTUBCS ) )**.8* 
$PRSTM**.3 

C COMPUTE THE LCG MEAN TFMPEPATURc C I F FERENC E FOR THE DE SU F ERH E AT S R 

LMTDCS=SC2/ ( HX FR CS* AP E ACS > 

C DETERMINE THE DESUPERHEATER METAL TEMPERATURE 

XXXXC2=EXP( ( SHOT-DSHOT l/LMTDQSl 
TM 50= (SHCT-XXXXC2*DSHCT ) / (1.0-XXXX0 2 1 
C DE SUP ERHFATER WATER SICE HEAT TRANSFER COEFFICIENT 

KW2=EC2/(TM50-TW1) 

C DESUPERHEATER STEAM SIC* HEAT TRANSFER COEFFICIENT 

KC2=EC2/(GD2**.8*LMTCCS! 

C SUPERHEATER STEAM SICE HEAT TRANSFER COEFFICIENT 

LMTD SF - ( TC2-TC1 ) /AUOG ( (TM20-TD1 )/(TM20-TD2) ) 

K C I = ER 2 / ( GO 1=** . £*LMTDSH) 

C CALCULATE THE TFPOTTLE VALVE FLCW CONSTANT 
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C 1 1 = GC4/ ( PD2*VENT) 

C SUPERHEATER OUTLET OE N S I T y 

CS*( FC1+P02 )/ (RC1+P.C2 ) 
C10=(FC2+PD3)/(RU2+RD3) 
RAYSFs(RCl+RQ2) /2.0 
RA VCS = ( RC2+RD3 ) /2.D 
C 7 = ( FC1-PD2)*RAVSH/GC1**2 
C9 S ( FC2-PD3 )*RAVCS/GQ2**2 



CIRCUATICN SYSTEM CALCULATIONS 



COMPUTE LOSS COEFFICIENTS FOR DCWNCCMERS 

FCP= ( l.C/d .74-2.C*ALCC10 (KSDDO) ) »**2 
C1CJ= <FCR*LCR/DCB+Z1+Z2+1 .0 ) / ( 2 . 0 * A C* * * 2*RW » 

FCS=(1.Q/(1 .74-2.0*ALCG19(KS0hC)) ) * * 2 
C 1 C 1 = ( FCS*LDS/00S+Z7+ZE+1 .C) / (2.0*ACS*+2*RW) 

ITERATIVE PR0C2CURE FCP BALANCING CIRCULATt.CN LCCFS 

SET INITIAL TR IAL VA LUES 

G C =2 C C . 0 
GS=20 C .0 
HH2=HH-EC2/G0 

START FIRST I TE RATI CN LOOP 

1 F H 6 = F H 2 +ED2 /GO 

X C = ( ER 2-GO* ( hFl-HW6) ) / ( G0*HFG ) 

FXF=FNCRM-HHTRCM 

HNBR=FXR*AMAX1 ( (HW1 -HW6 ),0.0)/<(HV,l + XC + FFC)-FV»6) 

L X R = F X R-HN BR 

R A 1= ( ( LXP./ ( X0*V FG ) »*ALCG((X0 4 VFG»/VF-*1.T )+RW*FNER )/HXR 
R C - 1 . C/( VF + XT'-* VFG) 

FR = ( 1 .0 / ( 1. 74-2 . 0*vALQG 10 ( K SDM3 I ) I **2 * < 1 . 0+1 6C78 .5 8* ( X0/ FC 1 » **.St ) 
. FSIMB* <FR*LR/CR + 1 .0 ) *FA1Z (2 .0 * AR * +2 *R C * *2 I +Z 3 / < 2 .0 *A R * * 2 *RW ) + 
JZ4/( 2.C*AR**2*RC) 

G0X=SCRT( (RW-RA 1 ) *HX R *G/ (C100+PSIMB ) ) 

400 FCR M AI (1H0»G15.5J 

IF CHECK VALUE EQUALS INITIAL GUESS, CONTINUE; ELSE REITEPATS 
IF (AES (GO-GOX ) .LT..01 )GO TO 2 
GC=( GCX-GC) /2 .C+GO 
GC TO 1 

START SECONC ITERATION LOCP 

2 CONTINUE 
XS s (CCl-GO*XC) /GS 
FXS = FNORM-FIFCR 

FNBS=FiXS*AMAX1 ( (HN1-HW2) ,0.0)/(<F>W1 + XS*FFG)-HW2) 

L XS=FXS-FN3 S 

RA2=( (LXS/(XS*VFG) )*ALCG( ( XS*VFG ) /VF+ 1 . 0 ) +RW*FNES ) / ( HXS ) 

R S 1 1 . C/(VF+XS*V C G) 

FS= ( 1 .0/ (1.74-2 .C*A LOG 10 ( KSOSC ) ))**2M1.0+16078.56MXS/FC1)**.S6) 
PS ISC = ( FS*LS/DS*1 .9)*PA2/(2.0*AS**2*RS*- < >2)+Z5/<2.C*4S**2*RW) + 
$Z6/( 2.C*AS**2*RS ) 

GSX=SCRT((Rli-RA2)*HXS*G/(C101+PSISC ) ) 

IF CHtCK VALUE EQUALS INITIAL GUESS, CONTINUE ELSE REITERATE. 
IF(ABS(GS-GSX) .LT..01ICO TO 3 
GS= ( GiX-GS ) /2.01GS 
GC TC 2 

START THIRO ITERATION LOOP 

hlft2X=(GW2/( GS+GC) )*(Hlr4-HV>l)+hWl 
C IF CHECK VALUE EQUALS INITIAL GUESS, CONTINUE ; E L £F REITERATE. 

I FCABS (HV»2-HV«2X ) .LT.. 1 IGO TO 4 
H K 2= ( FK2X-HH2 )/2 .0 + HW2 
GC TC 1 
4 CONTINUE 
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CALCULATE INITIAL CCNCITICNS 

G0P0=GC/PC 
GS&0=GS/RS 
MWIC=VCR*RW/2 .0 
H^IO^W 1*MWI0 
GXC0=hWl+X3*HFG/2.O 
CXS0=FWl+XS*hFG/2.0 
V R C J = (V0R/2.0) * R D 1 
HW60=HW 6 
GCI0=GC1 
G W R 0 = G 0 
G W 50 = G S 

RA IXC = (RAl*HXR-RW*HNiBR)/LXR 
RA2XQ= (RA2*FXS-RW*HNBS )/LXS 
LC1R=LAV0C 
L D 1 $ = L A VHO 
W P IT F ( 6 , 60 D ) 

60C F C FN A 1 (1HL) 

WRITE(6,INCCN1) 
wR ITS (6 ,601 ) 

6C1 FCF^Al C 1 H 0 ) 

WRITE(6,INC0N2) 

Vh F IT E 16,601 ) 

WRITE (6 , I NC CN3 ) 

WRITE (6,601 ) 

WRITE (fctCCWSTl ) 

WRITE (6 ,601 ) 

WP ITE (6 t C0NST2) 

UR IT r. (6 ,601 ) 

WRITE(6,C0NST3) 

WRITE (6,601 ) 

WRITE (6,C0NST4) 

WR IT E ( 6 , 60 1 ) 

WFITE (6 , CON ST 5) 

WPITE(6,6C1 ) 

WRITE ( 6 , CON ST 6 ) 

WRITE (6,603 ) 

WRITE(6,0LTPUT) 

WRITE (7 tINCO.Nl) 
WRITE(7,INCCN2) 

WRITE ( 7 1 1 NC CN 3 ) 

WRITE (7 , CONST 1 ) 

WRITE (7 , CCN £ T 2 ) 

WR ITE ( 7,CONST3) 

WFITE ( 7 , CCN ST 4 ) 

WRITE ( 7 ,CCiN$T5) 

W P IT E ( 7 , CONST 6 ) 
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